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Abstract 

We consider non-relativistic curved geometries and argue that the background 
structure should be generalized from that considered in previous works. In this ap¬ 
proach the derivative operator is defined by a Galilean spin connection valued in the 
Lie algebra of the Galilean group. This includes the usual spin connection plus an 
additional “boost connection” which parameterizes the freedom in the derivative op¬ 
erator not fixed by torsion or metric compatibility. As an example we write down 
the most general theory of dissipative fluids consistent with the second law in curved 
non-relativistic geometries and find significant differences in the allowed transport co¬ 
efficients from those found previously. Kubo formulas for all response coefficients are 
presented. Our approach also immediately generalizes to systems with independent 
mass and charge currents as would arise in multicomponent fluids. Along the way we 
also discuss how to write general locally Galilean invariant non-relativistic actions for 
multiple particle species at any order in derivatives. A detailed review of the geometry 
and its relation to non-relativistic limits may be found in a companion paper. 
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1 Introduction 

Newton-Cartan geometry was hrst developed shortly following the inception of general rel¬ 
ativity in an attempt to phrase non-relativistic physics in a manner that respects general 
coordinate invariance [1,2] and later developed further in [3-12]. Initial studies concerned 
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themselves with only the spacetime structure; however, recent applications in condensed 
matter physics have focused on covariantly coupling matter to this background. Perfect 
fluids in non-relativistic backgrounds were hrst studied in [13,14] while applications to the 
fractional quantum Hall effect have proven exceptionally rich [15-18], ranging from gravi¬ 
tational anomalies to energy transport. Newton-Cartan geometry has also naturally arisen 
in the study of non-relativistic holographic systems, where the boundary theory realizes a 
“twistless-torsionful” Newton-Cartan geometry [19-25]. 

One of the principle benehts of the geometric approach is that it furnishes a collection of 
background data consistent with Galilean invariance that may be perturbed to covariantly 
dehne currents and study linear response. These include a “clock” one form dehning a 
preferred notion of spatial vectors and elapsed time as well as a “spatial inverse metric” 
satisfying 


= 0 . ( 1 . 1 ) 

In the presence of a conserved particle current, one also has the option to couple to a 
background U{1) connection A^. In the case of a single Schrodinger held we then have 

S = J d^xVh Do'll; - ^Di'il;'^Dj'ijj^ ( 1 . 2 ) 

where — iA^, are the spatial components of the metric and we have taken 

= (^1 oj for simplicity. 

It was noted in [26] that for this action to be invariant under arbitrary space and time 
dependent changes of coordinates, the vector potential must carry the anomalous transfor¬ 
mation law 


5Ao = -i^dxAo - Axi^, 6A, = -^^dxA^ - Axdii^ - rnhij^ (1-3) 

where is an inhnitesimal diffeomorphism. Though simple in form, the hnal term is rather 
curious as A^ does not transform as a one form. The origin of this term was properly isolated 
in [27], identifying an implicit change in rest frame in the transformation above. The vector 
potential is then a true one form under diffeomorphisms, 

dA^ = -i^dxA^ - Axd^i^ (1.4) 


while under a Galilean boost in the manner 


(1.5) 
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where is a spatial vector representing the boost velocity. 

The backgronnd data is then the triple {n, h, A) np to a boost transformation given by 
(1.5). In this paper, we demonstrate that this is not the largest collection of data consistent 
with Galilean invariance and angment the backgronnd accordingly. Onr list of backgronnd 
helds inclndes 

( 1 - 6 ) 

Here is a veilbein transforming as a Galilean vector nnder local boosts and rotations and 
is eqnivalent in content to the clock form and spatial metric considered above, oj^b is a 
connection one-form valned in the Lie algebra of the Galilean gronp. It contains the nsnal 
spin connection 0 ;“^ = bnt also a “boost connection” = cj'^q which transforms nnder 
a choice of reference frame. We have also separated ont the single gange held appearing in 
previons treatments to two gange helds a and A so that we may describe independent mass 
(or particle nnmber) and charged cnrrents. It is then the mass gange held a that transforms 
in the manner (1.5), while the electromagnetic gange held A is boost invariant. 

This treatment also resolves a lingering qnestion of [27]. With the data (n, h, A) at hand, 
one may dehne a connection 

T^^i. = {df,Kp duhpp - dphp^) + n(^pF^)phP^, (1.7) 

where F = dA is the held strength of A. Here is a vector held parameterizing the 
choice of rest frame and is dehned by (2.22). This connection is boost invariant and so 
dehnes a sensible, frame independent geometry in the absence of torsion, bnt fails to do so 
generically. However, the Galilean spin-connection oj^b dehnes an invariant connection V 
on all backgronnds. In a companion paper [28], we consider the snite of possible constraints 
that may be placed on (1.6) in a boost invariant manner. In particnlar, the connection (1.7) 
may be obtained by the identihcation 


izr°' ACa = da (1.8) 

bnt we hnd that this is only possible when the torsion vanishes. 

One might reasonably hope that this issne may be avoided since torsion vanishes on a 
large class of physically relevant backgronnds^. However, stndying energy transport reqnires 
the introdnction of a Lnttinger potential $ [31], which arises in the spacetime approach as 
temporal torsion [18,32]. Thns if we hope to use Newton-Gartan geometry to study energy 
currents, we must know that our derivative operator respects Galilean invariance in such 

^Torsion can be relevant in the study of elasticity [29] and lattice defects [30]. 
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a background. To this end, we shall consider non-relativistic fluids coupled to (1.6) and 
perform an entropy current analysis and And results that do not agree with those of [33]. 

We begin in section 2 by reviewing the basics of Newton-Cartan geometry, boost transfor¬ 
mations and representations of the Galilean group that will be needed for our later analysis. 
Section 3 then introduces the mass gauge held a and boost connection m®. Before continuing 
on to fluids, we present in section 4 a brief detour through Galilean invariant actions and 
show how to write down actions for massive non-relativistic fields at any order in derivatives. 
The approach is seen to be equivalent to null reduction, but is intrinsic to the Newton-Gartan 
spacetime. 

Finally we develop first order fluid dynamics in our approach, beginning with a manifestly 
boost covariant presentation of the full set of diffeomorphism Ward identities in section 5. 
Section 6 then presents the fluid equations of motion and performs the entropy current anal¬ 
ysis. For multicomponent fluids, the results are summarized as follows (all coefficients are 
arbitrary functions of the temperature T, charge chemical potential /tq and mass chemical 
potential hm unless stated otherwise). The most general set of first order transport coeffi¬ 
cients in the parity even sector includes four sign semi-definite functions: a bulk viscosity, 
shear viscosity, conductivity and thermal conductivity 

((■ > 0, > 0, cte > 0, kt < 0- (1.9) 

as well as a thermo-electric coefficient 

ax. ( 1 . 10 ) 

The parity odd sector contains a Hall viscosity. Hall conductivity, thermal Hall conductivity, 
thermo-electric Hall coefficient, magnetization and energy magnetization 

r), dij, kr, ax, fh, fh^. (1.11) 

The magnetization determines the magnetic field induced pressure via the coefficient 

Jb = T‘^d,pdx + dgpdq + dppdM , ( 1 - 12 ) 

where p(e, q, p) is the pressure as a function of energy, charge and mass density. Kubo 
formulas for these coefficients are then presented. 
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In flat backgrounds with no Luttinger potential the constitutive relations are 

= p, p* = pw* f = q, + e, 

f = qu" + (Je{E" + Be'^Uj - Tc?Vg) + aEe"^{Ej + Bejku'" - TdjUq) 

+ ard^T + + e^WjTJi, 

^nc = + e + p - (9 - /bB^ m* - pa^'^Uj - pa^'^Uj 

+ T(Tt(^* + Be^^Uj - TaVg) - Tdre^^iEj + Be^ku’^ - TOjUg) 

+ ktO^T + kre^^djT - rhe^^Ej + BejkU^) + e^WjTh,, 

^nc = P^*“^ + (p - - fBB)g"^ - pa"^ - pa'-P (1.13) 

where is the mass current, the charge current, the energy current and the spatial 
stress and the fluid shear cr*-^ and expansion 6 are dehned by 

aij = diUj + diUi - Sij6, 6 = diu\ (1-14) 

The epsilon symbols are chosen with sign convention = ei 2 = 1. In summarizing these 
results we have made the choice of fluid frame (6.67). 

Note added: It was noted in [34] that our result misses a vorticity induced pressure 
coefficient due to a technical error. The correction may be found in (C.3) of their work. 

2 Galilean Symmetry 

The salient feature of non-relativistic physics is that of Galilean relativity, which, in its most 
familiar form, asserts that the laws of physics do not depend on a choice of inertial reference 
frame. These frames are related by the Galilean transformations 

t —)■ t x* —)■ Q'^jX^ — kH. (2.1) 

Here 0*^ G SO{d) determines the relative orientation and E G the relative velocity of 
inertial observers adopting coordinates (f, x*) and (P,x'®). 

In curved space one in general loses a notion of inertial frames and preferred coordinate 
systems. Rather, one can dehne only coframes 

ei, A = 0,l,---d ( 2 . 2 ) 

which form a local basis of 1-forms on the spacetime manifold Ai, which we take to be d -|- 1 
dimensional. The proper coordinate invariant statement of (2.1) is then that these coframes 
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transform as 



(2.3) 


Here spatial indices rnnning over the valnes 1,... ,d are denoted by lower case Latin letters 
a,b,... to distinguish them from spacetime indices A,B,.... Greek letters fi,u,... will 
represent coordinate indices while will be their spatial components. Of course, if 

the veilbein can be chosen to correspond to a global coordinate basis = dx^, we retrieve 
(2.1). Given a coframe we may of course dehne a set of frame vectors such that 

eX = = 6% (2.4) 

which transforms via the inverse of (2.3), 

(e„ e.) ^ (e„ e,) ■ (2.5) 

It’s easily seen from (2.3) and (2.5) that these transformations preserve the spacetime 
tensors 


n 





ft"" = 


( 2 . 6 ) 


The natural geometric setting for non-relativistic physics thus involves a positive semi- 
dehnite symmetric (2, 0) tensor and 1-form n^. These tensors are related insofar as 
spans the single degeneracy direction of 


= 0 . 


(2.7) 


These two helds are usually taken as the starting point in dehning a Newton-Gartan 
geometry, though we have chosen rather to go through a veilbein formalism since this will 
prove most convenient for our later analysis. They have clear physical interpretations. The 
“clock-form” dehnes a preferred notion of spatial direction at each point as well as an 
arrow of time: vector helds such that 


> 0 


being future directed. Any curve 7 also inherits a notion of elapsed time 


AT = 


n 


( 2 . 8 ) 


(2.9) 
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while serves as a spatial “inverse metric”^. 

The clock form defines a pointwise notion of spatial directions via vectors such that 
= 0. However, the Frobenius theorem tells us this notion may be integrated to a local 
codimension-1 hypersurface if and only if n A dn = 0. For our non-relativistic spacetime 
to carry an (at least local) notion of simultaneity, we must then demand that this holds 
everywhere. There are far more pressing reasons to take n A dn = 0 however, for a theorem 
due to Caratheodory [35] ensures the existence of closed timelike curves passing through any 
point violating this condition. To be precise, if n A dn 7 ^ 0 at the point p ^ AA, then there 
exists a neighborhood of p in which any two points may be connected by a future directed 
timelike curve^. We thus refer to spacetimes with n A dn = 0 everywhere as causal and will 
only consider such spaces throughout. 

In causal spacetimes, one may always choose coordinates such that n has no spatial 
components, and so the metric and clock form take the form 




( 2 . 10 ) 


where is everywhere a metric of Riemannian signature. 

2.1 The Galilean Group 

A few words on Galilean representations will prove helpful in what follows. The matri¬ 
ces appearing in (2.3) form a group under matrix multiplication called the Galilean group 
Gal{d). Throughout we shall refer to this defining representation as the vector represen¬ 
tation and denote matrices in this representation as Goframes then transform in the 

vector representation of Gal{d) and frames in the dual 



ca —t eB(A A- 


( 2 . 11 ) 


There is an equivalent definition of Gal{d) in terms of invariant tensor data. One may 


check that the matrices A"^^ are the unique ones that leave unchanged 




( 2 . 12 ) 


^The terminology can be deceptive. Since is degenerate it is neither invertible nor a metric. However, 
since it contains precisely enough data to define a unique Riemannian metric on spatial slices this terminology 
should not provoke undue confusion. 


^The cited theorem actually concerns null curves. However, it is easy to extend this result to timelike 
curves by adding a very small future directed component. 



This is the reason that a NC geometry contains precisely an invariant 1-form and degenerate 
spatial metric from a representation theoretic point of view, for the tensors 

= UAe^, (2.13) 

are then the primitive Galilean invariants that may be formed from the veilbein. 

Since Gal{d) is a subgroup of SL{d+ 1,M), the epsilon symbols eAo---Aa and with 

^oi - d = = 1 are also invariant tensors. We may use them to define a spacetime volume 

element 

(2.14) 

as well as a “raised volume element” 

^ gAo-AdgMo . . . gMd _ ^2.15) 

We caution the reader however that is not obtained from by raising indices 

since the metric is degenerate. These are both annihilated by the derivative operator 
that we shall define and in coordinates where (2.10) holds, 

Boi-d = V\h\e~^, 

We can also naturally contruct an invariant spatial volume element (with raised indices only) 


2.2 The Extended Representation 


Another representation of Gal{d) will prove remarkably useful when we turn to writing 
Galilean invariant actions (section 4) as well as in presenting manifestly covariant fluid 
dynamics. Explicitly, it is of the form 




j — 


/ 1 


0 o\ 

0G 0 


1 / 


(2.17) 


where = k°'ka = 5abk°‘k^- We shall refer to this representation as the extended represen¬ 
tation. It is d -|- 2 dimensional and has it’s own set of defining invariant tensors 


ni 


(l 0 o), 


jj 


/O 0 l\ 
0 5 “^ 0 

V 0 oj 


(2.18) 
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Here the indices I, J,... run from 0 to d + 1 and will always denote an object in the 
extended/dual-extended representation. The convention used in (2.17) is that the hrst 
row/column indicates the time component 0, the second indicates the full set of spatial 
components a and the hnal one the {d + 1)*'^ component. 

Although less familiar than the vector, the extended representation is often easier to work 
with since it admits the Lorentzian metric . Whenever working with this representation 
we will freely raise and lower indices with in the typical manner. Furthermore, any object 
in the extended representation may be projected to a Galilean vector via application of 

n^/ = o) (2.19) 

which is also invariant under the action of Gal{d). 


2.3 Timelike Vector Fields and Milne Symmetry 

We thus see that non-relativistic geometries are naturally dehned by a preferred clock-form 
and a metric that annihilates it. However, one thing that it cannot include is a preferred 
timelike vector held v^, for any such held is equivalent to a preferred notion of rest^. It is 
nonetheless often convenient to introduce such a for use in explicit formulas (for instance, 
in computing the Christohel symbols) and is essential for writing time derivatives (and thus 
any dynamical equation - see section 4). We will always normalize in the sense 


= 1 . 

Given such a we may then dehne a transverse projector 

P\ = < 5 ^ - 


( 2 . 20 ) 


( 2 . 21 ) 


whose upper index is orthogonal and whose lower index is orthogonal as well as a 
“lowered spatial metric” that “inverts” the degenerate to this projector 

Kuv'^ = 0, = P^. (2.22) 


In coordinates (2.10) we have 



(2.23) 


^Of course, this may be sensible in the presence of additional fields; for example, a background magnetic 
field or lattice establishes a preferred notion of rest with respect to which all velocities may be measured. 
However, the spacetime structure itself can make no such distinction. 
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Since has no physical signihcance, there must be a symmetry under shifts by a spatial 
vector so that (2.20) is preserved 


v'^ ^v>^ + F = 0 (2.24) 

This symmetry is often referred to in the literature as a “Milne boost” [36], and the projector 
and lowered metric also transform under this redehnition 

-)■ - k^Uu, -)■ - n^k^ - kf,n^ + k'^n^n^, (2.25) 

where = h^yk'^ and = h^yk^k^. Note that if n = e~^dt then the pullback of dehnes 
an invariant non-degenerate metric on constant t hypersurfaces. If we adopt a veilbein 
formalism for Newton-Cartan geometry, there is a natural timelike vector held in the form 
of Cq. Since Cq transforms under local Galilean boosts in the same manner as (2.24) it is 
natural to identify 


F = (2.26) 

and so Milne boosts and local Galilean boosts are also identihed. In the veilbein formalism 
we then also have 


P^u = h,, = (2.27) 

and the transformations (2.25) again follow precisely from (2.3) and (2.5). 

3 Non-relativistic Geometry 

To complete our description of non-relativistic spacetimes we also require a notion of parallel 
transport. This is specihed via a connection V that preserves the above data 

= 0, VaF" = 0. (3.1) 

In the Lorentzian case, metric compatibility completely determines the connection (up to 
torsion). Unfortunately, as is well known, (3.1) does not do so in a Newton-Gartan geometry 
and there are many distinct connections (even if we require vanishing torsion) all satisfying 
these conditions. 

Various authors have dealt with this ambiguity in different ways. In the recent condensed 
matter literature, this extra freedom has been hxed by introducing additional conditions to 
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further constrain the connection. In [16] this was done by demanding the curl-freeness of 
the vector held 

h\[fj,Vu]V^ = 0. (3.2) 

The associated Christoffels are 

{df,Kp + d^h^p - dphp ^), (3.3) 

where we have assumed torsion-freeness for simplicity. Of course, since and hp^ are not 
a boost invariant quantities, the connection so dehned is not either. 

As mentioned in the discussion below (1.7), recent works [27] have taken the approach 
that this freedom should be hxed by the held strength of the gauge held coupling to 
particle current, asserting rather 

h\[pS^ u]V^ = --Fpy. (3.4) 

giving Christohels 

= v^d(^pny) + {dphyp + dyhpp - dphpy) + (3.5) 

where one of the Fpy indices has been raised with a h^''. With the transformation law 

Afj. Ap + kp — —Upk'^, (3-6) 

this dehnes a boost invariant connection in torsion-free backgrounds. 

However, as discussed in [27] and reviewed in section 1, this is lost upon the introduction 
of torsion. The resulting notion of parallel transport then depends on a notion of rest frame 
and so is inconsistent with the principle of Galilean relativity. In [28] we consider the possible 
restrictions consistent with boost invariance that one may place on the connection and hnd 
limited freedom when the background is arbitrary. We thus take the point of view that the 
connection should not be restricted when studying energy transport. 

3.1 Newton-Cartan Geometry 

In our approach then, a Newton-Cartan geometry is a d -|- 1-dimensionaI manifold At with 
a collection (n, h, V) of a 1-form a rank-d positive semi-dehnite metric h^'^, and a con¬ 
nection V such that 

hP^riy = 0, WpTiy = 0, = 0. (3.7) 
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The connection is then simply extra data above and beyond what appears in the clock-form 
and spatial metric and is boost invariant by dehnition. 

The additional data may be parameterized in a boost non-invariant way by introducing 
a v^. It’s derivative is some tensor 


= A/ (3.8) 

that includes information on the acceleration, shear, expansion and twist of the vector held 
v^. The Christoffel symbols are then 

{d^hyp d^h^p — dphpy) -\- — (T^ 

(3.9) 

where we have also included possible non-zero torsion These Christoffel symbols were 

recently derived in the formalism of Koszul connections in [37]. For a veilbein approach 
see [28]. Here VLpy = involves the vorticity and acceleration of which contains the 

data unhxed by the torsion and (3.7). We shall refer to this as the Newton-Coriolis 2-form. 
One may check that under boosts, O transforms as 

+ 29[^(/Ci,] - ^n^]k‘^) - {kx - ^nxk‘^)T^pu (3.10) 

so that (3.9) is invariant even in the presence of torsion. In deriving (3.10) we have made 
liberal use of the identity Vxhpy = —2AA(^nj,) which is an immediate consequence of the 
dehnitions. 

In relativistic physics we work with Lorentzian geometry where the torsion may be chosen 
independently of the metric, but it is well known that the analogous statement is not true 
for non-relativistic physics. Anti-symmetrizing the expression VpU^ = 0 gives the condition 

nxT^^ly = {dn)py. (3.11) 

Although we expect dn = 0 on most physically relevant non-relativistic spacetimes^, it is 
useful to keep this torsion around at intermediate stages when studying energy transport. 
For this reason it is important to know that our geometry is sensible in the presence of 
torsion. 

3.2 A Galilean Spin Connection 

The data contained in a Newton-Cartan geometry can be equivalently stated purely in a 
veilbein formalism. In this approach, rather than the clock form Up and spatial metric 

^dn = 0 is equivalent to a notion of absolute time. For details we refer to the discussion in [28]. 
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we have the coframe transforming as a Galilean vector. The complete data of the 
connection is then equivalent to a Galilean spin connection oj^b dehned by 

^ = CbUJ^^a (3.12) 

which implies the usual transformation law for connections 

^d{^ + A"^C'd(A (3.13) 

for an arbitrary local Galilean transformation A^b- The curvature and torsion are then 
dehned in the usual way 


R'^b = doj^B + A Bi = De^ = de^ + oj^b A (3-14) 

and are boost covariant. 

The compatibility conditions (3.1) then imply 

0 = e\V^el = 0 = e“e^VA(e^e^h''^) = 20;^^“'^- (3.15) 

That is, they are equivalent to oj^b being in the Lie algebra of the Galilean group gal{d) 


Ul B — 


0 0 


ZU UJ i, 


(3.16) 


where ujfj_ab = ^fi[ab]- The spin part of the Galilean connection is familiar. It dehnes the 
connection on a spatial slice and may be used to covariantize actions involving helds that 
transform under local rotations through the covariant derivative 

i/j ^ e = dfj, + —ujfZ^Jab- (3-17) 

where Jab sue the generators of rotations. 

The boost part ci7“ is, however, new. It transforms as a vector under rotations but as a 
connection under boosts 


zu^ ^ zu^ + Dk‘^. (3.18) 

where Dk°^ = dk^" + u°'bk^- It is the boost connection that parameterizes the freedom in the 
metric compatibility conditions and is equivalent to the tensor A^^" considered earlier 

A/ = = (3-19) 

It is then our point of view that a non-relativistic geometry is specihed by the pair uj'^b), 
which is equivalent to the treatment of the previous section. 
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We shall occasionally find it usefnl to consider the spin-connection valued in the extended 
representation 


/ 0 0 o\ 

j = I w°' uj°‘b 0 
\ 0 -Wb Oy 


(3.20) 


Note that the I, J indices here indicate a homomorphism of from the dehning to the 
extended representation of the Lie algebra, and not application of any of the invariant tensors 
considered previously. 


4 Non-Relativistic Actions 

We now have the necessary tools to present non-relativistic actions in a manifestly invariant 
manner. We begin with the Schrodinger action in d dimensional flat space 

s= y d'‘*'x f , (4.1) 

where — iqA^ is the electromagnetic gauge-covariant derivative. The obvious 

diffeomorphism covariant generalization of this expression to curved space is 

s = y j (4.2) 

where the integral includes an implicit factor of the volume element e. For a spinful 
Schrodinger held, the covariant derivative will of course include the spin connection as men¬ 
tioned above. One immediately retrieves (4.1) by going to hat space 

= (l o) . fr = . (4.3) 

and selecting a frame where =^10^ . 

Unfortunately, this action makes explicit reference to and so is not manifestly boost 
invariant. This problem necessarily arises whenever time derivatives are involved since the 
only diheomorphism covariant expression that includes a cIq is However, the action 

can be made to be invariant by introducing a gauge held for mass, 

Df, = dfx- ima^ - iqA^ (4.4) 
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that transforms in the manner 


a^, + kf, - -n^k'^. (4.5) 

As presented this may seem ad hoc, however, a direct non-relativistic limit naturally 
identihes with the lapse and shift and so corresponds to gravitational forces (see [28]). The 
transformation (4.5) then naturally follows from an ambiguity in the ADM decomposition 
of a Lorentzian spacetime. In [27], this gauge held was identihed with the electromagnetic 
held. Ehectively, this amounts to collecting a and A into the combination a + ^A, and 
so the treatment could only describe particles whose mass was tied to their charge. By 
separating from we have two independent symmetries 17(1 )m and U{1)q and have 
the added beneht that we may describe non-relativistic systems whose components have 
multiple charge to mass ratios. We will exploit this added freedom in section 6 to generalize 
the analysis of [27] to huids that have independent mass and charge currents. 

The action (4.2) is boost invariant and a clear generalization of the hat space action. 
However it is a rather curious combination of one derivative terms, two derivative terms, 
and frame dependent quantities with precise factors of m to make everything work. It 
would be nice to know how to systematically generate such combinations at any order in 
derivatives. To address this problem, we begin by rephrasing (4.2) in a manifestly invariant 
manner. The result (4.9) is essentially equivalent to a null reduction and may also be found 
in equation (3.23) of [27]. However, this will teach us how to select the proper Galilean 
invariant combinations of time and space derivatives at any order we choose. 

We shall eventually wish to work with multiple particle species, so let M denote the 
generator of U{1)m and Q denote that of U{1)q, which we will take to be anti-hermitian. 
For our single held we have for instance M^jj = imip and The derivative 

D^ = df,- a^M - A^Q (4.6) 

is U{1 )mxU{1)q gauge covariant, but not boost covariant and the time and space derivatives 
Dq = Da = e^D^ transform as 

Do ^ Do + k^Da - Da ^ Da- kaM. (4.7) 

These can however be collected into a derivative valued in the dual of the extended repre¬ 
sentation (2.17) 

Dj = (^Dq Da , so that Dj —)■ (A“^)'^/Dj. (4.8) 
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(4.9) 


The obvious two derivative quadratic action is then 

S = J Di'ip^D^'ip 

which upon expanding yields precisely (4.2). The form clearly mirrors that of null compacti- 
hcation, where plays the role of the higher dimensional metric and (4.8) is the derivative 
operator in that space. The identihcation D^+i = im follows by retaining only modes of 
momentum m in the compactihed circle. The advantage of this presentation however is it 
compactly describes the action (4.2) in a manner intrinsic to the non-relativistic spacetime. 

This is a simple rewriting; the true beneht of collecting and Da'ip into Diijj is that 
it’s now clear how to generalize the precise cancellations in (4.2) to arbitrary boost invariant 
actions of any order in helds and derivatives. To do so, we extend the dehnition (4.8) to a 
derivative operator Vj on tensors in the extended representation 

Vju'^ = Diu'^ + uj/ ku^ , T^iVj = DjVj — vk(^i^ j, ( 4 - 10 ) 

and the obvious generalization for more general extended-valued tensors. By the hrst term 
we mean simply apply the operator (4.8) to vj. 

In the dehnition (4.10) we have converted the spacetime index on to an extended 
index via the projector (2.19) 

Uk^ j = j- (4-11) 

One may then check by hand that Vj is locally covariant, which follows in the usual way 
upon realization that j = j, giving the required cancellation. For instance 

'Di'u'^ = Dju'^ -\- ojj'^ 

(A jDk^A"^lu^) + (A mu^ + (A Ld^{A ^A^mu^ 

(A jA ^+ (A jA'^ lojk^mu^ 

+ {A-Yi {Dk{A\)u^ - n^Ke^^d^{A\)u^) 

= {A-YiA\Vku^. (4.12) 

It’s now a simple exercise in index contraction to write down Galilean invariant actions 
of any order k in derivatives. Given a collection of helds in some representation of 
U{1)m X U{1)q (to account for dihering charge/mass ratios), they are of the form 

F = S[^l^\ ..., P,, • ■ ■ Vi4, (4.13) 

where all indices have been contracted using the invariant tensors and nj and of course 
all terms have U{1)m x U{1)q charge zero. 
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As an example, consider the term 




(4.14) 


In flat space this gives rise to 


/ -imDQijA 


= Tr 


/ imD^xlj imDQ'ip\ 


imD^'ip D'^D^ip D^Dq'iJj 
yirnD^i/j DoD^i/j DoDoi/j J 


DaD.iP^ -zmDaip^ 
y—irnDoip^ —imDcip^ 

= m^^-DQDo'ip^ + 2Doi>^Doip - DoDopj) 

+ im{{Do, DaPip) + DaD^iP^D^D^P j 


(4.15) 


which one may explicitly check to be invariant under global boosts. The spin connection is 
not necessary to form globally invariant flat space actions, but is needed to get the proper 
local boost invariant action, supplying the necessary terms involving derivatives of via 

(3.19). 


5 Ward Identities 

In preparation for the fluid analysis of section 6, we now turn to the Ward identities asso¬ 
ciated with diffeomorphism invariance. These merely express local momentum and energy 
conservation and serve as the equations of motion for thermodynamic variables upon taking 
the hydrodynamic limit [38]. So far as we can tell, they were first derived in this language 
in the work of [39], which specialized to flat geometries and latter augmented to include a 
curved metric hij in [40]. 

The full diffeomorphism covariant Ward identities on arbitrary backgrounds were first 
spelled out in [32], but the treatment was not boost invariant, including a stress tensor and 
energy current that depended on a choice of v^. This was resolved in [27], which first defined 
a boost invariant, diffeomorphism covariant stress-momentum tensor and stated its Ward 
identity without reference to v^. The work-energy equation however can only be stated in a 
boost invariant way in the presence of a boost invariant vector held for instance in the 
presence of a fluid [33]. 

Our approach is essentially equivalent other than issues concerning the connection already 
discussed, but we take the time to restate these identities in a language that will prove useful 
to us later on. The only mild innovation is that we collect the stress tensor, energy current, 
momentum current, and mass current into a single object carrying an index valued in 
the extended representation of the Galilean group. The dependent Ward identities then 
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manifest themselves as a single invariant equation (5.8) for work done by external fields on 
T^I- 


5.1 The Stress-Energy Tensor 

To do this we begin by comparing (2.17) and (4.5), noting that the veilbein and mass 
gauge held may be collected into a single object valued in the extended representation 



which we shall call the “extended veilbein”. As discussed for the derivative operator in 
section 4, one can think of the Newton-Cartan geometry as a null reduction of a Lorentzian 
geometry in one higher dimension [27,28] and the “extended veilbein” is simply the reduction 
of the higher dimensional veilbein. Care should be taken with this terminology since there is 
no sense in which this forms a veilbein on spacetime as the collection is necessarily linearly 
dependent. 

The other background helds are then the Galilean spin connection uj^b and the vector 
potential A. We dehne their associated currents as 

5S = j • (5.2) 


is then the charge current and the spin current. The tensor carries an index 
valued in the extended representation and transforms covariantly. We can then easily re¬ 
trieve the dependent treatments of previous works by parsing out this tensor into it’s 
components. 

The currents as dehned previously in the literature [32] 

(jS* = y + f5A^ (5.3) 


rely essentially on a choice of but now are seen to naturally arise in the single object 





-Pncb 

__ a I 

nc ^ r nc y 


(5.4) 


Here and in what follows the label “nc” will be affixed to those currents measured in the “lab 
frame” dehned by v^. is the response to the variation of the purely spatial components 
of the veilbein and so corresponds to the spatial stress. Pnc6 is the momentum current, is 
purely spatial and as we shall see is equal to the spatial mass current pncb when matter is 
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not charged under boosts, a fact that we will assume in the remainder of this section. We 
shall refer to as the stress-energy tensor for convenience, despite the fact that it contains 
data beyond energy and stress. 


5.2 Diffeomorphisms 


In the presence of massive charged helds, the background helds are then an extended veilbein 
e^, a spin connection oo^b and the electromagnetic gauge held A. Their transformations 
under inhnitesimal diffeomorphisms ^ are the standard ones 

6e^ = L^de^ -|- = Dc^e^ -|- , 

doj^B — i^doo^B T di^uj^B — I'^R^b T Di^uj^Bi 

SA — L^dA “t" di^A = -t- dt^A^ ^5.5) 


where i^F represents the contraction etc. We have dehned an extended torsion tensor 


= De^ = de^ -|- j A e'^ 



(5.6) 


includes the full information on the spacetime torsion dehned in (3.14) as well as a 
“mass” torsion / = da — Wa /\ e°‘. We shall see that / couples to the mass current as an 
external held strength and so should not be thought of as a Newtonian gravitational force, 
which is encoded within the covariant derivative (see [28] for a more detailed discussion of 
Newtonian gravity in relation to our geometry). 

The transformations (5.5) diher from covariant expressions by a U{1)q gauge transfor¬ 
mation —l^A and an inhnitesimal local Galilean transformation l^oj^b- Since our theory is 
assumed to be gauge and local Galilean invariant we may perform a simultaneous diheomor- 
phism, U{1)q and Gal{d) transformation, under which we have 


6e^ = Di^e^ + 5uj^B = i^iR^B 5A = i^F. (5.7) 


Variation of the action then yields the Ward identity 

-e^iD, - T\,)r", + = F^,f + R^a,.s’^b^. (5.8) 

As ever, represents the locally Galilean covariant derivative D^r^j = 

One might worry about the appearance of in the hrst term. However, if mass is 
conserved (as we shall always assume), the explicit term drops out, as it multiplies the 
Ward identity 


(V^-TV)P^ = 0, 


(5.9) 
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which follows in the usual way upon a U{1)m gauge variation 5a = (IKm- Of course, charge 
conservation follows as well from 5A = dAg 

(V^ - = 0. (5.10) 

Despite the deceptive appearance of the mass gauge field, we shall find (5.8) very useful 
owing to it’s manifest boost invariance. 

Equation (5.8) contains both the work-energy equation and momentum conservation 
equation of previous works. To get a sense of this equation, the reader may hnd (6.31) 
helpful, where these components are isolated and rendered in a more familiar form for fluid 
systems. Finally, we note that since we have defined the stress and spin currents in a veilbein 
formalism, our stress is a response to a shear at hxed which necessarily introduces 

torsion = De^. The physical stress tensor rather measure response to a shear at fixed 
spatial torsion. The proper redefinition of currents will introduce additional terms into the 
Ward identity involving double derivatives of the spin current and so will not enter into our 
calculations with the power counting scheme we will adopt. For this reason we shall also 
drop the final term in the Ward identity. 

5.3 Local Boosts 

It is well known that Galilean symmetry fixes the momentum to the flow of matter. The 
relevant Ward identity to demonstrate this is that of infinitesimal local boosts, under which 
we have 


5a = kac'^, = -Pe°, 5w^ = Dk^ = dk^ + uj\k\ (5.11) 

while the spin connection is neutral. Variation of the action then gives 

0 = j {p^ka - Pl^ka + s\^Dk^) , (5.12) 

from which we have 

pSc = <0 - (D, - (5.13) 

Interestingly we see that the oft-cited non-relativistic identity between the mass and 
momentum currents may be violated in the presence of matter that couples to the boost 
part of the spin connection 

Df, = d^- af,M - A^Q - (5.14) 
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Here Jab are the generators of rotations and Ka of boosts, which, together with translations 
Pa, time translations H and mass M, close to form the Bargmann algebra [41] 

[jafe^ jcdj ^ ■ (^^acjbd _ ^adjbc _ ^bcjad ^ ^bdjac^^ ^ 

p^] = i - S^^P^) , [J^\ K^] = i , 

[P^,K^] =[H,K^] = -iP^. (5.15) 

The identity = p“(, is then generically violated whenever there is matter that trans¬ 
forms nnder the sub-algebra spanned by {Jab, Ka} in which Ka is represented non-trivially 

^ (5.16) 

Such representations were hrst considered in the work of Levy-Leblond [42] and enumerated 
up to and including spin 1 in [43,44]. However, since we do not know of any condensed 
matter systems that realize these representations we shall assume Ka = 0 in what follows so 
that = 0 and the identity between momentum and mass currents is retained 

pSc = P™- (5.17) 

6 Non-Relativistic Fluids 

In this section we turn to non-relativistic fluid dynamics to illustrate the advantages of our 
approach. The program of fluid dynamics is to describe the fluctuations of thermodynamic 
variables in slightly out of equilibrium media. These variables include the temperature and 
velocity as well as a chemical potential for each conserved charge. For us, there are two 
such quantities, the electric charge and mass, whose associated chemical potentials we will 
denote /ig and /im (/^m being the same boost invariant hm dehned in [33]). In all we have 
d -|- 3 degrees of freedom 

T, idQ, Pm, u^, (6.1) 

where the velocity has been normalized so that = 1. The equilibrium properties of 

the system are then completely characterized by an equation of state such as p(T, fig, ^m), 
which expresses the pressure as a function of the state variables. The entropy, charge, mass, 
and energy densities are then dehned by 

dp = sdT + qdfiQ + pdpM, e + p = Ts + qpg ppM- (6.2) 
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Just out of equilibrium, the thermodynamic variables are allowed to vary slowly in space 
and time, where slow is compared to the mean free path and mean free time so that a local 
equilibrium is always a good approximation. The Ward identities 

(V^ - T\,)f = 0, (V^ - T%,)ff = 0, 

- el{D, - (6.3) 

are then sufficient to serve as equations of motion since they are the same in number as the 
thermodynamic variables. We need only constitutive relations for the currents in terms of 
these degrees of freedom and their derivatives. 

Since we are perturbing away from equilibrium, these constitutive relations naturally 
organize themselves in a gradient expansion where higher derivatives of (6.1) take on di¬ 
minishing importance. To complete our power counting scheme one needs also specify the 
backgrounds to be perturbed around. In this paper we shall assume a trivial background in 
equilibrium so that 

V^, (6.4) 

are all hrst order in derivatives. This is not a unique selection and corresponds to a choice 
of regime in which we expect our results to be applicable. One may for example consider 
backgrounds of large constant magnetic held in which B appears at zeroth order and the 
analysis will be signihcantly altered. 

The constitutive relations so obtained are not in general sensible and could lead to an on- 
shell decrease in entropy, the canonical example being that of a negative shear viscosity [38]. 
As such we must also impose the second law of thermodynamics as an additional constraint 
on huid hows, diminishing the freedom present in the gradient expansion and returning a 
reduced set of transport coefficients. This procedure has been carried out in many systems 
including 2 + 1 dimensional relativistic normal huids [45] and 3-1-1 dimensional relativistic 
superhuids [46]. 

In this paper we will work only to hrst order in derivatives to illustrate our method in the 
non-relativistic setting. It’s well known that even at hrst order, the dynamics is very rich 
in the parity odd sector when d = 2 and so we shall eventually restrict to two dimensional 
huids. 

Although our calculation is very much along the lines of [33], there are two key diherences. 
The hrst is the existence of an independent mass current, a necessary element for treating 
systems with multiple components. One might for example consider mixtures of He-3 and 
He-4 in their normal phases. It shouldn’t be surprising that such systems admit a richer 
transport sector. In this work we hnd that a multicomponent huid admits one additional sign 
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semi-definite transport coefficient, one additional unconstrained parity even coefficient and 
three additional parity odd coefficients compared to the single component fluid. A succinct 
overview of these results may be found in section 6.8. Our second point of divergence is that 
stressed in section 1 and we will be using an connection without the kinematical restrictions 
imposed in previous works. 


6.1 Covariant Currents 


As we’ve already noted, the currents come naturally assembled into a single stress-energy 
tensor r^/. However, in the presence of a boost invariant fluid velocity we may go further 
and define Galilean frame invariant notions of energy and stress separately. Heuristically, 
this corresponds to defining them to be as measured in a frame co-moving with u^. In 
this section, we demonstrate the details of how to do this. Our approach in this regard is 
essentially equivalent to that found in [33], though we restate it here in our language. These 
currents at hand, we parse the Ward identity (5.8) into something a bit more familiar: the 
work-energy and Navier-Stokes equations. 

To begin, consider the fluid velocity as measured with respect to some lab frame 

= ( ^\ 

The existence of a preferred allows us to define P^b = — u^ub and a h^B such that 

h^'^hcB = P^B unambiguously®. We shall also make use of 


(6.5) 


P^i 


P^bT^^i 


/ 0 0 o\ 

5\ 0 ) 


( 6 . 6 ) 


which projects extended indices to transverse vector indices. 

Any vector index be decomposed uniquely into a part parallel to and perpendicu¬ 
lar to Ha, while an dual vector index may be decomposed into a part parallel to ua and 
perpendicular to u^. For instance, for a vector and covector wa we may write 


= vu^ + 


Wa = wua + w'a 


(6.7) 


where v'^ua = = 0. 

®In this section only, these tensors correspond to the fluid frame, not the frame defined by v^. 
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We should like to perform a similar decomposition for extended indices I for which we 
have the preferred vector 


= g^-^nj = 


/0\ 

0 

VV 


( 6 . 8 ) 


However, is a null vector of the Lorenztian metric gu and is thus perpendicular to itself, 
making the above procedure impossible. In Lorentzian geometry one usually continues by 
introducing a second null vector whose inner product with is 1, though such a f is 
not unique. One can then decompose any index into parts parallel to parallel to , and 
perpendicular to both. Luckily in the presence of a background fluid, there is a preferred 
way to select such a vector. We will dehne in the extended representation so that is both 
null = 0 and projects to the fluid velocity in the sense . In components, 

the so dehned is of the form 


u = 


/ 1 \ 


(6.9) 


and automatically satishes nju^ = 1. 

Using this extended velocity vector, we now decompose the lower index of into parts 
parallel to ri/, parallel to uj and perpendicular to both 


- p\i + 


( 6 . 10 ) 


where by dehnition = t^ju^ = 0. When an /, J,... index is transverse in this sense, it 

may always be written as the pnllback of a nniqne tensor with a transverse A, B,... indices. 
In our case we have t^i = where bu^ = 0. We now continue, decomposing the 

upper index of t^B into a part parallel to and a part perpendicular to ua 


t B = u^Pb + 


B- 


Altogether we have 


T^i = e^rii - p^ui - {u^Pb + T^b)'A^ 


I- 


( 6 , 11 ) 


( 6 . 12 ) 


We can invert this dehnition to hnd 

p^ = -T^in\ / = -ubt^iP^^, (6.13) 
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These are the boost invariant currents we shall use in the fluid analysis and can be so deflned 
whenever there exists a preferred velocity vector to draw on. They correspond to the energy, 
mass current, momentum and stress as measured by an observer comoving with the fluid. 
The relation to the component decomposition (5.4) is 

^nc - + |Pnc“7 ’ \PIc) ’ 

rjnAB ^ ^ 

yO ~ p^U^'v}’ 

where we have used the Ward identity We note that is simply the Milne 

covariant energy current deflned in [33]. 

We now turn to restating the Ward identity (5.8) in terms of these currents. To simplify 
matters, we shall now take so that 0^ and define the transverse projector 

and inverse metric accordingly. 

A few words on torsion are due before we proceed. Throughout we shall restrict to 
backgrounds in which there is no spatial torsion. Spatial torsion is necessary in the study 
of media with dislocation defects, where the presence of torsion simply indicates a nonzero 
Burgers vector [30,47] (for a classic discussion of dislocations in elastic media see [29]) and 
would be an interesting element to include in future work. For our purposes however, the 
fluid analysis is greatly simplified by discarding it. 

The minimal temporal torsion (3.11) may also be set to zero by assuming the clock-form 
is closed. Although this is the case in all physically relevant situations, temporal torsion 
is useful for studying energy transport [48] as it is equivalent to coupling the system to 
a Luttinger potential [31]. As noted previously, our formalism also allows for a non-zero 
“mass” torsion f = da — A which will be zero on physical backgrounds. Although 
including mass torsion in our analysis will not prove as fruitful as temporal torsion, there is 
no essential difficulty in doing so and so we keep it around. 

The boost covariant manner in which we set spatial torsion to zero but keep these effects 
is a decomposition of the extended torsion tensor along the lines of the previous section, 
keeping only components parallel to and 

= u^G + n^j (6.15) 

where G = dn and 7 = / + UaT°‘ — ^v?G. This decomposition would in general obtain a 
third term that is both n and u orthogonal corresponding to spatial torsion but which we 
have set to zero here. The spacetime torsion tensor is then = u^G^y, whose pullback to 




- P “ 


= P^bP^. 


(6.14) 
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a slice is zero if we assume nAdn = 0. The Ward identity for diffeomorphisms then becomes 


-el{D, - G,)t^i = (6.16) 

We see that G serves as a held strength coupled to energy and 7 a held strength coupled to 
mass. We have also dehned G^ = —G^^u’^ = 

Now supplement the decomposition (6.12) of the stress-energy with one of the extended 
veilbein 

el = n^c^ + u^b^ +q^n^ +^ (6.17) 

where by dehnition 

Uiq^ = njq^ = 0 , UiQ^ ^ = njq^ ^ = 0 , = 0 - (6.18) 

A straightforward computation then shows that 

“I . (6,19) 

The term drops out of the equation of motion as it multiplies the continuity equation 
for mass. The rest then reads 

- + q^^)iD^ - G^)t’'i = + F^yf - G^yS’'. ( 6 . 20 ) 

Contracting with we have 

u\D^ - G,)t^i = e,ff + E^f + G^e^ 

{V, - G^)e^ = + E,j^ + G,e^ + r^iDyU^. (6.21) 

where 

= 'Jf.uu'', = Ff.yu’'. ( 6 . 22 ) 

are the held strengths observed by co-moving observers. The hnal term is rather mysterious 
looking but can be easily evaluated 

r'^iDyU^ = {s’^rij - p’^uj - (pEpA + T"'A)'^'^i)DyU^ 

= -{u'^PA + T^A)^^lDyU^ 

= -{u^PA + T’^A)DyU^ 

= -(uVa + T^)V,u" 

= -Pua'^ - - \0g^^Tyx. (6.23) 
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In the second line we have used that Dyirii = 0, nju^ = 1 and ujDyU^ = 0 and in the third 
line that = 0. The acceleration, shear and expansion of the fluid appearing in the 

above formula are dehned as 

a^^ = 9 = (6-24) 

Altogether 

(V^ - G,)e^ = (e^ - + G,e>^ - (6.25) 

which is the work-energy equation, including work done by the external helds as well as 
dissipated by fluid shears. This matches exactly the covariant work-energy equation as it 
appears in [32] and [27] except for the additional mass current that couples to 7 ^ 1 ,. 

To obtain the Navier-Stokes equation we raise the /i index on (6.20) 

-(V, - G.)(r^g^^) = YuP'^ + (6.26) 

The stress tensor term on the LHS is evaluated as 

= {u^pA + T^A)n^/g"M = wV;. + T%. (6.27) 

Similar manipulations to (6.23) gives 

= p^UiD^q^^ + {u'^px + T^)V.(e^n^g^'^). (6.28) 

Now e\ll^iq^^ = so that the hnal term drops out. To evaluate the hrst term we have 


UjD^q^f^ = Ui{V,yq^>^ + ojJ jq'^^) = 

/ 0 0 0 \ / 0 \ 


= 0+001 






= -eM = -V.nA 


(6.29) 


VO/ 


\ 0 -w^b 0/ 

Rearranging the vector equation of motion and using the mass conservation Ward identity 
we have 


(V, - G,){u^P\p^ + - p^V,u^ 

(V. - G,){p^u'^ + mV" - pu^u" + = l\p" + - G\e'^ (6.30) 

where p = n\p^. 

In summary, we may restate the Ward identities in terms of the covariant currents as a 
work-energy equation and Navier-Stokes equation, plus conservation laws 

(V^ - G,)e^ = (e^ - a^)p^ + + G,e^ - 

(V, - G,)(p'^m" + mV" - pu^y" + T^") = 7Vp" + ^Vj" - 
(V^ - G,)ff = 0, (V^ - G,)3^ = 0. 
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(6.31) 






6.2 Perfect Fluids 


The Ward identities (6.31) serve as dynamical eqnations once constitntive relations have 
been snpplied, specifying the cnrrents and in terms of the thermodynamic 

degrees of freedom 

T, /ig, Pm, u^- (6.32) 

At zeroth order in derivatives, the most general tensors we can constrnct nsing these variables 
and the Newton-Cartan strnctnre are 

j^ = qu^, p^ = pu^, = (6.33) 

Here g, p, e and p are fnnctions of (T, pg, Pm)- They are identihed with the thermodynamic 
charge density, mass density, energy density and pressnre and so satisfy the relations (6.2). 

Feeding these into the Ward identities (6.31) we obtain the perfect flnid eqnations of 
motion 


q + q9 = 0, p + pP = 0, e+{e + p)9 = 0, 

pa^ = pe^ + qE^ + (e + p)G'^ - V>. (6.34) 

In these eqnations and those that follows dotted objects indicate the material derivative, 
/ = The hnal eqnation is simply Newton’s second law and an obvions covariant 

generalization of Enler’s eqnation [38]. It expresses the fact that flnid particles will tend flow 
along geodesics, deviating only dne to the exertion forces from external backgronnd helds 
and internal pressnre. 

6.3 Fluid Frames 

To go beyond perfect flnids one needs to expand the cnrrents to hrst order in derivatives. 
The constitntive relations are then the perfect flnid ones pins 0{d^) corrections 

f = {q + Q)u^^ + i^^, p^^ = {p + q)u^ + 

= {e + 8)u^^ + e, = {p + V)h>^'' + 71^'^. (6.35) 

In the above the vector corrections are dehned to be transverse 

= n^p'" = = 0 (6.36) 

while the tensor correction is traceless 

(6.37) 
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(tt^^ is of course already transverse since T'^'' is). This is convenient as it separates the hrst 
order corrections into irreducible representations of SO{d). 

The decomposition (6.35) is subject to a well known ambiguity stemming from the need 
to dehne T, pg, and out of equilibrium. Any such dehnition is admissible so long 
as it reduces to the equilibrium values at order zero and so is subject to a d + 3 parameter 
0{d^) held redehnition 

T —)■ T + 6T, fiQ —)■ fiQ + SfiQ, fiM hM + d/iM) —>■ + 6u^ (6.38) 

called a huid frame transformation (not to be confused with a Galilean frame transforma¬ 
tion). To deal with this ambiguity we may either £x the frame by imposing extra conditions, 
or work in a manifestly frame invariant manner. 

Frame transformations are worked out in [33] (see [49] for a relativistic treatment) and 
in this section and the next, we refer the reader to this treatment for the details. For our 
purposes we only note that (besides those related to the entropy) the complete set of hrst 
order frame invariants is 

S = V - d^pS - dqpQ - dppg, 

= (6.39) 

P P 

Although we shall usually take p to be a function of temperature and the chemical potentials, 
here we have taken p = p(e, q, p) and the partial derivatives c?e, dq and dp are defined accord¬ 
ingly. Note we have an additional vector frame invariant compared to either the relativistic 
case or non-relativistic single-component huids due to the presence of the conserved current 

p^. 

6.4 The Entropy Current 

It is convenient to separate out the hrst order entropy current into a “canonical part” and 
corrections. Here the canonical part is dehned to be that combination of currents chosen to 
match the equilibrium identity (6.2) 

= pu^ + - Pgf - PmP^ 

i'S. = + + + (6.40) 

Out of equilibrium the entropy how will in general deviate from the canonical part 

+ C- (6.41) 
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This separation is helpful since is a frame invariant, and so is as well. It’s divergence 
is a quadratic form in first order data 

(V„ - 1(e„ - TV, (^)) - T£»‘(v„r - tg,). 

(6.42) 

6.5 Constitutive Relations 

All the necessary tools are now available to carry out the analysis outlined at the beginning 
of this section; first write out the most general constitutive relations for the first order 
frame invariants and then impose the second law of thermodynamics. We shall take d = 2 
throughout. The formulae will prove somewhat simpler if we instead take our independent 
variables to be 


T, Vq—^q/T, 

In terms of uq and um, the thermodynamic identities read 

dp=^-^dT+ Tqdi'Q + Tpdi'M, = s + I'Qq + pmP- 

Now consider the available first order data 


V^T, V/,z/Q, 






G 


flip • 


Separating into irreducible representations of S'0(2) we have 


where 


Scalar 

Vector 


Data 

e b B u 
(T) (f'g) {j>m) 

(«^) 


Symmetric Traceless Tensor 


a 


flU 


= rX 

= u'^VX- 


B = 
b = 
e = 


(6.43) 


(6.44) 


(6.45) 


(6.46) 
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Recall that by we mean the “spatial volume element” which is boost invariant. 

Not all this data is independent on-shell. We may thus use the Navier-Stokes equation to 
eliminate one vector degree of freedom and one scalar each for mass conservation, charge 
conservation and the work-energy equation. The eliminated data is indicated by parentheses 
in the above table. 

The above amounts to a decomposition of of the form 

= n^a'' + (6.47) 

Had we used the restricted connection such that 2ga[i_iVu]V°‘ = —F^y not all this data would 
be independent, but would be identihed with the electric held and uj with the magnetic 
held. 

Also note that we have not included a “torsional magnetic held” G = in the list 

above since this is zero on causal backgrounds. One could certainly include this and compute 
away, but we do not particularly trust our results when there is no notion of time evolution. 
In particular the second law of thermodynamics would be essentially meaningless. In either 
case, the point is moot as including a nonzero G does not introduce further constraints. 

Finally, the most general hrst-order constitutive relations for the frame invariants con¬ 
sistent with spacetime symmetries are 

S = -C9- hb - JbB - flu, 

+ geE^ + ggG^ + utV^T + gqV^vq + 

+ GeC^ + GeE^ + GgG^ + GtV^T GqV^Vq + 

KeE^ -\- KqG^ -\- 1 

+ + KeE^ KgG^ + -|- Kq'V^I'q -f 

= {Ce9 + lb + + ^^u)u^ 

+ Cee^ + CeE^ + CgG^ + CrV^T + CqV'^z/q + CmV'^z/m, 

+ l~e^ + CeE>^ + CgG^ + CrV'^T + CqV^z/q + CmV'^z/m, 

(6.48) 

In this we have dehned the “dual” operation 

= £^^''vy, (6.49) 

on vectors and symmetric two tensors. It has the properties 

(6.50) 

We have similarly used tildes to label parity odd response coefficients, e.g. ge is the normal 
electrical conductivity and ge is the Hall conductivity. 
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6.6 Entropy Current Analysis 

Now let’s move on to determining those constraints that result from imposing the second 
law of thermodynamics 


(V^ - G^)s^ > 0. (6.51) 

Start by considering the genuine second order data in the entropy production 

(V^ - G^)Cl2_a = + Cbh + (Cb - Ce)B + 

+ CrV'T + CqVVq + CmWm + CeV^e'^ + > 0 

where we have used the identities 

- G,)E, = -B- B9, = 0. (6.52) 

The hrst is just Faraday’s law and is equivalent to the closedness of E (note we do not 
have a corresponding identity for 7 ). The second similarly follows from dG = 0 as well as 
n A dn = 0. 

We conclude that all the coefficients listed above must vanish and so only has contri¬ 
butions from the remaining parity odd part 

= CbBu^ + CbE^ + CgG^ + CtV^T + CqV^z/q + (6.53) 

The rest of the divergence of is then 

(V;, - G^)C = - (rdepdrCB + ^dgpdQCB + ^OppOmCb'^ BO 

+ OtCbE^V^T + OqCbE^V^vq + OmCbE^V^vm 
+ (Cr + ^tCg)C'^V^T -|- (Cq -|- OqCg)^^'^+ (Cm + 

+ {BqCt - dTCQ)V^TV^UQ + {SmCt - dTCM)V^TV^UM 
+ (^mCq ~ ^qCm)V'^i^qV^i^m (6.54) 
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which is supplemented by the canonical entropy production 

+ ^aE{E, - TV^uq){E^ - TWuq) - ^kt{V,T - TG"^)(V^T - TG^ 

+ T^fhbO + —fsBO + —fijUjO + —aeC^E^ + —KeC^G^ + — (ug + Ke) E^Gfj, 

- - aee^V ^,vq + ^ E^V + ^aME^V ^um 

+ — (Jg^ G^V^uq + —i^mG^V^vm 

- (^T + V^TV— ^KmV'^TV^Z/m — O'm'^^VqS/ fj,VM 

+ - ^(ug - 

- - deC^V^^UQ -^(^aT + 

- ^de + —dg^ E^V^uq — —^mE^'^ fj,i'M 

~ y ~ ^uq — —HmG^'^^vm 

- ^dT - V^TV^Vq + ^KmV'^TV^Z/m + dMV'^Z/gV^l/M (6.55) 

In the above we have made the identihcations —y:crg = (Tg and —|:Kg = so that the 
entropy production due to electrical and thermal conductivity factors into a perfect square^. 

^In fact, what one should do is demand that the quadratic form defined by these transport coefficients 
be degenerate and positive semi-dehnite (degenerate so that that equilibrium solutions exist in non-zero 
background fields). This immediately gives these identities. 
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Demanding the second law then requires 


C > 0, ?? > 0, ue > 0, kt < 0, dc = Kg, 

CTg d? 

L = fb = 0, fs = T'^depdrCs + dqPdqCB + OppOmCb, 

\ 

Cq + dgCc 


b-e — —o'g — ~-^bq — Tax, 


\ / ^kx + i^kq^ 


^QtCb ^ 


^ + ^Be^ 

= kE + ^kq 

1 

dqCB 

= 

de + ^dg 

/ V / 


^QmCb } 


^ M J 


^dgCu — OmCo^ 



SmCt — dr Cm 

= 


\ dxCq — dqCx ) 


\^kq - ax) 


To untangle the differential constraints, begin by dehning 

f = Cb, = Ct + i^tCg, Thq = (q + dgCc + T f, ThM = CM + dMCG- 


(6.56) 


(6.57) 


These then read 


^ fix ^ 


^ rp2 “t” -pS ^ 


(dxf\ 


^ rp(TjC' ~\~ Jp2^E^ 

hq 

= 

^Re + ^Rq + / 

1 

dqf 

= 

O'E + 

\hM ) 


^ J 


\dMf) 


y M J 


/ 


^dju — OmJ 
1 


T_ 

h^hM — 


\ 


rp I J'S 

y T^ hQ -\- -^kq H“ dx f ^ f J 


^dqllM — dMhq^ 

Om^t — dxhM 

\ dxhQ — dgllT ) 

This leads to several consistency relations on the four functions (6.57) 

/ dohM — (0 \ 


'qbm 
dMhx — dxhM 

\ dxhQ — dqhT J 


0 

\V 


(6.58) 


(6.59) 


The hrst comes from comparing the hnal component of the second equation to the hrst 
component of the third while the second follows from comparing the hnal component of 
the hrst equation to the second of the third. The hnal condition results from combining 
the second component of the hrst equation, the hrst component of the second and the hnal 
component of the third. The vector hq Jim^ is then curl free as so must be the gradient 
of some function g{T, uq, em) 

fix = dxg, hq = dqg, Iim = duQ- (6.60) 
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6.7 Summary of Results 

This solves the full set of restrictions imposed by the second law. Before summarizing results, 
the following redefinition of transport coefficients will simplify the final answer 

T/^m, 

ar ^ d'T + Orfh, Rt ^ + OTm^. 

(6.61) 

after which frame invariants are 

5 = -(0 - fsB 

= ge - TV^vq) + - TG^) + ge {e^ - + GriV^T - TG>^) 

— rhG^ + V^m 

S^^ = Tgt {E>^ - TV^^uq) + - TG>^) - Tgt {e^ - + Kr(V^T - TG^) 

— rhE^ — 2m^G^ + ^^171^. (6.62) 

The most general set of first order transport coefficients are then as follows. There are 
four sign semi-definite functions of all three thermodynamic variables 

C > 0, > 0, o'E > 0, kt < 0- (6.63) 

These are the bulk viscosity, shear viscosity, conductivity and thermal conductivity, all of 
which are zero for dissipationless fluids. One sign indefinite parity-even coefficient exists, a 
thermo-electric coefficient 


gt (6.64) 

which determines the charge flow due to thermal gradients and the energy flow due to 
electromagnetic fields. 

The parity odd sector is richer, including six unconstrained parity odd coefficients 

V, &E, Rt, dr, rh, rhe, (6.65) 

a Hall viscosity. Hall conductivity, thermal Hall conductivity, thermo-electric Hall coeffi¬ 
cient, magnetization and energy magnetization, rh and the equation of state determine the 
magnetic field induced pressure 

Jb = T^d.pdT + dqpdQ + dppdM ■ (6.66) 
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The attentive reader may note that the parity odd response to thermal gradients differs 
from the parity odd response to the Lnttinger potential by the energy magnetization (and 
similarly for the thermoelectric Hall coefficient), as compared to the well-known resnlts 
of [31,50]. This is because these works assumed vanishing equilibrium currents, as pointed 
out in the footnote below equation (4.10) of [50]. In general, pure curl persistent equilibrium 
currents may arise, given by the magnetizations, and in this case the proper relationship is 
that given above (6.62). 

To get a feel for these results, it’s helpful to £x a fluid frame and write the results for the 
non-covariant currents dehned in (5.4). We choose our frame so the physical mass, charge 
and energy correspond with the thermodynamic ones and the velocity is that of the mass 
current 


g = ^ = ^ = 0, = 0. (6.67) 

The frame invariants and are then simply the hrst-order deviations and Using 
(6.14) to retrieve the non-covariant currents from this data we have 

= p, p" = pu" f = q, £nc = + e, 

f = qu'' -f aE{e^E^ + Be'^^Uj — TcIVq) -|- aE£^^{e'^Ej + BsjkU^ — TdjVQ) 

+ aTe^<9*(e-'^T) + (e'^T) + e'^£*^(e-'^m), 

4 c = + e + p-C 0 - IbB^ 

+ Tar(e^U* + Be^^Uj - T^Vq) - TaTe'^e^Ej + Bsjku’^ - TdjVq) 

-h KTe^d''{e~^T) + Rxe-^dj{e~^T) — me^^{e‘^Ej + BejkU^) + e^'^e*'^5j(e“^*rTie), 

4ic = pu^u^ + (p ~ ~ — fja^^ (6.68) 

In (6.68) above and in what follows U® is dehned to be the electric held in the lab frame 
so that the comoving electric held used earlier is® 

0 -E,\ / -E,u^ 

EiBeijjyu^J \e'^Ei +BeijU^ 

9 and cr®-^ are the curved space quantities dehned in (6.47) . Using the connection (3.9) we 
®We hope the careful reader will forgive the notational dissonance. 


(6.69) 
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find that 


+ P\u^ = + P\u^ + e‘^F*o, 

= - h!^ + 

^ = e'^V*(e-V) + e^V^'(e-‘^MO “ (6.70) 

where V* is the covariant derivative on a spatial slice and 

6 = — 7 ^ — dJVhe~^u^) = e^'Vi(e~^u^) + -e^0hij. (6-71) 

Vhe-'^ ' ' 2 " 

6.8 Results for a Single Component Fluid 

A single component fluid satishes additional constraints since the charge and mass currents 
are proportional and we investigate these constraints in this section. Since this is also the 
case of applicability for [33] it will allow for a direct comparison of our results. 

Let the single constituent be of charge e and mass m. The charge density and mass 
density are then related to a single function n, the number density 

q = en, p = mn (6.72) 

and the thermodynamic relation (6.44) takes the form 

dp = + Tndh> (6.73) 

where v = evq + mvM is the total chemical potential. All thermodynamic functions must 
be a function of this combination. 

Now since mv^ = the vector frame invariant is zero which gives the restrictions 

cr^; = (Tt = = 0, m = fh(T), dr = —rh', Tdx = —rh. (6.74) 

The latter two relations imply that ??f is a linear function of T and so that dy is a constant 
independent of the thermodynamic state variables. 

This simplihes to 

= kt{V>^T - TG>^) + kr(V'^T - TG^) + (6.75) 

Since we cannot have dependence on vq outside of the combination z/, we must also have 
d'p = 0. 
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In summary, for a single component fluid there are three sign semi-dehnite functions 
including a thermal conductivity and two viscosities 

C > 0, h ^ 0) kt < 0 (6.76) 

which exhausts the parity even sector. The parity odd sector contains three free functions 
including a Hall viscosity, thermal Hall conductivity and energy magnetization 

?7, Rt, fhe- (6.77) 

In terms of these coefficients, the frame invariants are 

5 = -Ce 

= KTiV^T - TG^^) + RTiy^T - TGf^) - + V^m,. (6.78) 

One may equivalently state this result in terms of a conductivity cte and a Hall conductivity 
Re by exchanging the hrst order data — TG^ for eE^ — TV^Ty — ma^ using the equations 
of motion (6.34). 

Finally, the noncovariant currents in the frame (6.67) are 

f = q, f = qu\ ^nc = + O 

+ e + p- (9^ E - 

+ KTe'^9*(e“'*’T) + RTe^e'‘^dj{e~^T) + 

T^i = pEu^ + (p — C9)h^^ — ~ rja^^. (6.79) 

This differs from the results (1.13-18) of [33], which includes two parity odd parameters 
besides the Hall viscosity and thermal Hall conductivity and a differential relationship de¬ 
termining the magnetic held induced pressure, a transverse Ei term and a curl diT in terms 
of them. 

We have also checked the results for lowest Landau level huids [48] with the alterations 
to the derivative operator discussed here and found that they survive without change. This 
can be understood in part in terms of our equations of motion, which differ from [27] as the 
acceleration is independent data and not tied to E^. Projection to the lowest Landau 
level (which may be thought of as a massless limit) removes the acceleration terms in Navier- 
Stokes equation (6.34) and we reproduce the constraint equation 

V> = +(e + p)G^ (6.80) 
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used in [48]. The only change to the calculation would then be in the available hrst order 
data, which is augmented in our approach. However, a detailed calculation does show that 
the new terms drop out after the entropy current analysis and we retrieve the previous 
results. 


6.9 Kubo formulas 


The transport coefficients we have found are quite familiar and have been subjected to 
extensive study in the literature and calculated for a number of systems. Calculation from a 
microscopic theory usually proceeds by the use of Kubo formulas. The techniques to derive 
these are now standard and Kubo formulas for all the transport coefficients presented above 
have been given in the literature. Here we present them in our notation and in the frame 
(6.67) for the readers convenience. 

The relevant retarted correlators, including contact terms are 

([r«W,T“( 0 )]>, 

Gy(x) = . 

Gs-u) = (^) - "(^”) , 

G%"[X) = - iOixO) (|t>‘(x),t>'(0 )]) . (6.81) 

In terms of these, the viscosities are 


C = — lim 




2iuj4 


rj = — lim 

on-O 


(<^+) 


2iuj4 


= _ lim (6,82) 

2iuj^ 


where we have introduced the projectors 

Y^ijkl _ n{kH)j _ }ixijxkl 
2 


^ (6.83) 


and u:+ = uj + i5 for a small, positive 5. We recommend [51] for a careful computation of 
these Kubo formulas. 

The equations for the conductivity and thermoelectric conductivities are prototypical 
examples and hrst found in [50,52] 




aE = lim 

u ;—>-0 Z % U ). 




Oe = lim 

cj —>-0 




Tar = lim 

2 icj. 




Tar + fh = lim 

1 *1—>-0 2iuj. 


(6.84) 
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Kubo formulas for the thermal conductivities were hrst computed in the classic work [31] 
where the Luttinger potential was introduced. We hnd them to be 


Tkx = 


lim 

cj—^O 


(‘^+) 

2ia;+ 


TRt + 2me 


lim 

oj—>0 


2iu+ 


(6.85) 


As discussed in section 6.7 the themoelectric and thermal Hall conductivities differ from the 
parity odd response to the chemical and Luttinger potentials by magnetizations m and 
respectively, unlike the standard formulas found in [31,50]. The Kubo formulas for ax and 
Rx are therefore completed by expressions for the magnetizations, derived in [48]: 


m — T dxfn = — lim 

|kK0 




2me — T dxThe = — lim 

|kK 0 




(k) 


( 6 . 86 ) 


7 Outlook 


In this paper we have considered the most general geometric backgrounds consistent with 
local Galilean invariance and developed the theory of hrst order dissipative huids on such a 
background. This formalism at hand, there are a number of prospects for further investiga¬ 
tion. One direct application would be to perform the huid analysis carried out here for more 
general systems. Non-relativistic superhuids and superhuid/normal huid mixtures as arise 
for instance in partially condensed superhuid Helium and it would be interesting to see the 
restrictions imposed by Galilean symmetry. 

We have also presented a program for writing down invariant actions of massive helds 
that realize the Galilean group linearly. These are not however the most general actions 
consistent with non-relativistic symmetries. For instance, when matter is charged under 
boosts V' one may write down a non-relativistic form of the Dirac equation that 

is linear in both time and space derivatives [44]. 

One might also consider systems with spontaneously broken symmetries in which the 
Galilean group is realized non-linearly on a collection of Goldstones. It’s long been under¬ 
stood how to write down the most general actions for a non-relativistic Goldstones to lowest 
order in derivatives [53] but we believe use of the extended representation will prove useful 
in constructing actions to any order. 

It would also be instructive to consider the most general ehective actions one may write 
in terms of the background helds o, A). Gonsider for example the (universal sector 

of the) ehective held theory of the fractional quantum Hall ehect [16,17,54], 


S = 


1 

47r 


[vA A dA + Ku f\ dA -|- k'o; A doo) 


(7.1) 
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where oj = \e°‘’^ujah- This is a perfectly sensible effective action for the background helds 
(A,e“) and useful for studying electric and viscous transport. 

However, studying massive transport also requires coupling to the U{1)m gauge held 
a. Since the microscopic action for a minimally coupled single component system always 
contains A and a in the combination A + (see equation (4.6)) one might expect that the 
correct effective action contains only this function of A and a. The resulting action is however 
not boost invariant and so physically unacceptable. Upon identifying boost transformations 
with the anomalous diffeomorphisms of [26], this problem is essentially the one considered 
in [15] where the problem is solve by improving the action order by order in derivatives so 
as to impose invariance. It would be useful to have a manifestly geometrical way to write 
this term and so generate these corrections using the technology developed here. 
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